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We present a phenomenological theory for spatiotemporal chaos (STC) in Rayleigh-Benard con¬ 
vection, based on the generalized Swift-Hohenberg model. We apply a random phase approximation 
to STC and conjecture a scaling form for the structure factor S{k) with respect to the correlation 
length ^ 2 - We hence obtain analytical results for the time-averaged convective current J and the 
time-averaged vorticity current fl. We also define power-law behaviors such as J ~ ~ and 

^2 ~ where e is the control parameter. We find from our theory that p, = 1, v > 1/2 and 
\ = 2p p V for phase turbulence and that p = 1, v > 1/2 and X = 2p + 2v for spiral-defect chaos. 

These predictions, together with the scaling conjecture for S{k), are confirmed by our numerical 
results. Finally we suggest that Porod’s law, S{k) ~ l/f, 2 k^ for large k, might be valid in STC. 

PACS numbers: 47.54.-fr, 47.20.Lz, 47.20.Bp, 47.27.Te 


I. INTRODUCTION 

Rayleigh-Benard convection (RBC) has long been a paradigm in the study of pattern formation . This system 
consists of a thin horizontal layer of fluid heated from below. There are three dimensionless parameters to describe 
the system |^. The Rayleigh number R = ga(PAT/Kv is the control parameter, in which g is the gravitational 
acceleration, d the layer thickness, AT the temperature gradient across the layer, a the thermal expansion coefficient, 
K the thermal diffusivity and v the kinematic viscosity. Under the Boussinesq approximation, only the density of the 
fluid is temperature dependent. Then the Prandtl number <j = vj k is all one needs to specify the fluid properties. 
The third parameter is the aspect ratio T = L/d where L is the horizontal size of the system. When the Rayleigh 
number R surpasses a critical value Rc, the fluid bifurcates from a static conductive state to a convective state, in 
which the velocity profile u = (uj_,Uj;) and the temperature-deviation profile 9 form certain self-organized patterns. 
Those patterns also depend on the boundary conditions (b.c.) at the horizontal surfaces of the container. The two 
most studied b.c. in the literature are the rigid-rigid b.c., under which the fluid cannot slip, and the free-free b.c., 
under which the fluid does not experience any stress. 

The patterns and the corresponding stability domain in RBC have been studied extensively in the classical work 
of Busse and his coauthors [^j|] in the (i?, a, k) space with k the wavenumber. This stability domain has hence 
been known as the “Busse balloon” in the literature. For rigid-rigid b.c., the parallel roll state is predicted to be 
stable inside the Busse balloon for any given a. Surprisingly, more recent experiments and numerical studies 
|p|-|l^ revealed that for a ~ 0(1) and T > 50 the parallel roll state yields to a spatiotemporally chaotic state even 
for states inside the Busse balloon. This spatiotemporally chaotic state, called spiral-defect-chaos (SDC), exhibits 
very complicated dynamics both temporally and spatially . Its discovery has since stimulated many experimental 
[P,[7 | JTi|Ji^ , theoretical |Q and numerical efforts to understand it. Despite these efforts, few insights have 

been obtained so far. 

For free-free boundaries, Zippelius and Siggia and Busse and Bolton found that the parallel roll state 
is unstable against the skewed-varicose instability immediately above onset if cr < 0.543. Busse et al. ||l^ further 
investigated the dynamics involved and conjectured a direct transition from conduction to spatiotemporal chaos (STC). 
This spatiotemporally chaotic state is called phase turbulence (PT). Recently we reported a large scale (F = 60) 
numerical simulation of the three dimensional hydrodynamical equations for a = 0.5 under the free-free b.c. H- 
From that simulation, we confirmed the direct transition to PT above onset and studied various properties of it. The 
patterns we found have very complicated spatial and temporal dependences. 

The dynamics of the two types of STC in RBC (SDC and PT) is very complex. From the theoretical point of 
view, it is far from clear whether methods developed in studying ordered states, such as the Galerkin method |^||, 
the amplitude equations [T^ , the phase dynamics [^ , etc., can be helpful at all. It seems that new concepts and 
new theoretical tools are needed in studying STC ||20|| . We know that thermodynamic systems can be characterized 
to a certain degree by global quantities such as temperature, pressure, density, etc.. A natural question one may ask 
is: Is it possible to characterize STC by some global quantities? The answer to this question is unknown at present 
but evidence for a positive one is very encouraging. Recently two of us found from extensive numerical studies that 
global quantities such as the time-averaged vorticity current and the spectrum entropy have different behaviors in 
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the parallel roll and SDC states |^. Furthermore, we demonstrated in our recent paper on PT that both the 
instantaneous and the time-averaged behaviors of global quantities carry valuable information about the state. Even 
if global quantities may not be sufficient to describe all the dynamical details of STC, as a first step, we believe that 
knowledge of such quantities will broaden our understandiM of STC. 

The generalized Swift-Hohenberg (GSH) model of RBC [£l|-g3| is widely accepted for theoretical study. This model 
is derived from the three-dimensional hydrodynamic equations, but is much simpler to study both numerically and 
analytically. After the corresponding z-dependences are separated, the vertical velocity Uz{v, z, t) and the temperature- 
deviation field 9{y, z, t) are reduced to an order parameter ^/’(r, t) in two-dimensional space r = (a;, y), while the vertical 
vorticity • V x uj_(r, z, t) is reduced to u;z{r, t) pl|-p^. There are two coupled equations in the GSH model, one for 
'0(r, t) and the other for a;z(r, t). The convective patterns in RBG are completely determined by the order parameter 


V'(r, t). Numerical solutions of this model or its modified versions have not only reproduced most patterns observed in 
experiments but also resembled experimental results relatively well P|, |T^ , P|J^ ]. But there are some shortcomings 
in the model ||^,^: The stability boundary of the model does not coincide with that of hydrodynamics; it induces 
an unphysical, short-ranged cross roll instability; and both the shape and the peak position of the power spectrum 
for SDC are different from those in the real system. Even so, owing to its simplicity and its qualitative resemblance 
to real systems, this model is very valuable in studying RBC. 

In this paper we present our analytical calculations, using the GSH model, of the time-averaged convective current 
J = A~^ f dr'i/’^(r,t) and the time-averaged vorticity current ft = A~^ f drujl{r,t) in both PT and SDC, where F{t) 
represents the time-average of F{t) and A is the area of the system. We carry out our calculations in Fourier space 
so the total number of modes '0(k, t) considered is infinite. By assuming the time-averaged two-point correlation 
function C'(ri,r 2 ) = '0(ri, t)^/>(r 2 , t)/V'^(r, t) is translation invariant in STC, i.e., C'(ri,r 2 ) = C{ri — r 2 ), we find that 
the phases of two fields are uncorrelated in time unless they have the same wavenumber k. Furthermore, 

we apply a random phase approximation (RPA) to STC in which four-point correlation functions are approximated 
by products of two-point correlation functions. Using this RPA, we derive J and D in terms of the time-averaged 
and azimuthally averaged structure factor S{k) = ^*(k, t)'0(k, t)/J. We further assume that S{k) obeys a scaling 
form kS{k) = ^ 2 tF[{k — kmax)^ 2 \, in which ^2 is the two-point correlation length, F{x) is the scaling function and 
kmax is the peak position of kS{k). Applying this assumption, we obtain explicit formulas for both J and Q. More 
precisely, we find that J = Joe— which depends on unknown but experimentally measurable parameters, where 

e = {R — Rc)/Rc is the reduced control parameter and Jo and J{ are both known. On the other hand, we find that 
H = wi J^/^l for rigid-rigid b.c. and D = ui 2 wJ ‘^/^2 for free-free b.c., where w is related to the width of the scaling 
function F{x) and is experimentally measurable. The other coefficients uji and 102 are known exactly. Furthermore, 
by assuming power law behaviors such that J ~ e^, ^2 ~ and D ~ e^, we predict from our theory that /i = 1, 
V >1/2 and \ = 2pL + 2v for rigid-rigid b.c. and ^ = 1, ^ > 1/2 and \ = 2pL + v for free-free b.c.. This prediction and 
the scaling assumption for S{k) have been verified by our numerical solutions for both PT and SDC. 

Our paper is organized as follows. In Sec. we introduce the GSH model in Fourier space and derive the basic 
formulas governing the time-averaged convective current J and the time-averaged vorticity current Q for any pattern 
in RBC. We present a simple theory for PT in Sec. HI. This theory is based on a conjecture of Busse et al. 0] 
that PT can be described by an infinite number of modes lying on a single ring in fc-space with different orientations. 
Although the result for J seems to agree with our numerical result extremely well, this is somewhat accidental since 
it predicts, incorrectly, U = 0 and also neglects the strong couplings between modes of different fc’s. In Sec. IV we 
introduce the RPA of our theory and use it to calculate explicitly J and U for both PT and SDC. The results are 
expressed in terms of the structure factor S{k). Sec. ^ includes three parts. We first postulate the scaling form of 
the structure factor S{k) and expand both J and ft in the leading order of 1/^2- We then define power law behaviors 
for J, ^2 and U in PT and compare the results from our theory and our numerical work 0, which agree very well 
for both the exponents and the amplitudes. In the last part of this section we define power law behaviors for J, ^2 
and D in SDC and test our theoretical formulas with our numerical results ||^ ]. The agreement between our theory 
and our numerical work is good in general, except for the amplitude of D. In Sec. ^we discuss the large-A: behavior 
of the structure factor S{k). We conjecture that Porod’s law in which S(F) ~ for large k with a 

characteristic length, might be valid for STC. In the last section, we summarize our results and discuss some related 
issues. 


II. BASIC FORMULAS 

Near the conduction to convection onset, the velocity field u(r, z, t) = (u^, Uj,) and the temperature-deviation field 
0{r,z,t) in RBC can be approximated by pl|-|^ 
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where V is the gradient operator in two-dimensional space r = {x,y). For both free-free boundaries at z = 0,1 
and rigid-rigid boundaries at z = ±i, the explicit forms of uo{z), wq{z) and 9o{z) are given in Ref. one takes 
Co ( 2 :) = 1 for free-free boundaries and Co ( 2 ) = z"^ — \ for rigid-rigid boundaries. Notice that the vertical vorticity 
• V X Ux is now replaced by Co(- 2 )o 22 (r, t) in which ^^(r, t) = — V^C(i’) ^)- Inserting Eq. (|l|) into the three-dimensional 
hydrodynamical equations in RBC and applying a few more approximations, one ends up with the two-dimensional 
generalized Swift-Hohenberg (GSH) model of RBC |2l| - |^ . Although some of the approximations are not systematic, 
the amplitude equations for the GSH model and the hydrodynamical equations are the same in the leading order near 
onset. There are two coupled equations in the GSH model, one for the order parameter '0(r,t) and the other for the 
mean-flow field C(r, t). The convective patterns are completely determined by the order parameter The GSH 

model has been proven very successful in characterizing convective patterns under quite broad conditions . 

In the GSH model, the order parameter t) satisfies 

TO [dti^ + U • V^] = [e - {ej^kl){V^ + klf] ^ - N[^]. (2) 

Here iV['0] is the nonlinear term to be specified soon and U(r) is the mean flow velocity given by U(r) = VC(r, t) x e^, 
in which 


[dt - - c^)] V^C = • [V(VV) x . (3) 

In the GSH equations, the reduced Rayleigh number e = (R/Rc) — 1 is the control parameter, in which R and Rc 
are the Rayleigh number and its critical value at onset. The Prandtl number a parameterizes the fluid. While 
is the critical wavenumber at onset, the other parameters model the prope rties of the system. The values of these 
parameters depend on the boundary conditions (b.c.), more precisely l22y29|] , 

for free — free b.c. : Rc = 277r'^/4, kc = TrjV^, tq = 2(1 -|- cr“^)/ 37 r^, 

= 8/3^2, = 6 , = 0; (4) 

for rigid - rigid b.c. : Rc = 1707.762, kc = 3.117, tq = (1 -F 0.5117ct-^)/19.65, 

= 0.148, gm = 24.77, c" = 10. (5) 

It is easier to analyze the GSH equations theoretically in Fourier space than in real space. By convention, we define 
the Fourier transformation and its inverse transformation of an arbitrary function F'(r) as 

F{k) = jf dre-**^ "F(r) and F(r) = ^ F(k)e*‘^ ", ( 6 ) 

where A is the area of the system. Note that F*(k) = F'(—k) for any real function F{r). It is easy to check that Eq. 
can be rewritten in Fourier space as 

Todtipik) -h E(k) = r(e; k)'tp(k) - N{k), (7) 


where V (r) = tqU • Wip and 


Since is real, one has '(/j*(k, t) = ■)/((—k,t). The nonlinear fV(k) term has been evaluated at onset | |2^ , 

^(k) = X! ■ k2)'tp*{k2)ip{k3)-ip{k + k2- ka), (9) 

k2,k3 

where the coupling constant 5 (cos a) is given in Ref. with a the angle between k and k 2 . Rigorously speaking, 
the exact forms of Eqs. (|)> i) and (^ are derived near onset and deviations from them in real physical systems are 
possible for large enough e. But we disregard such complexity and take them as our model for further study. 

One may take an adiabatic approximation {dt = 0) in Eq. (H) by neglecting the first term on the left-hand side. 
This term is small in comparison with the other terms, which can be verified by applying the same perturbation as 
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that in phase dynamics H- With this approximation, it now is easy to solve Eq. (||) for C(k), which indicates that 
the mean-flow field is slaved by the field. We are also interested in the vertical vorticity ^^(r) = — V^C(r). 

From Eq. (||), it is straightforward to get that 

d)^(k) = /(k; k2)-0(k2)Vi(k - ka), (10) 

k 2 


where, with an exchange of index k 2 ^ k — k 2 , 

/(k; k 2 ) = g 2 ) “ 2k • k 2 )(e^ • k 2 x k). (11) 

Applying these results, one may easily evaluate the mean-flow contribution to Eq. (Q), which is given by 

^(k) = Y ^^(k;k2;k3)'0*(k2)'!^(k3)V)(k-l-k2 - ka), (12) 

k2,k3 


where 


u(k; k2;k3) 


gmTo [e^ ■ k X (k 3 - k 2 )][e^ ■ ka x k 2 ](fc| - fc|) 
2 cr |k3-k2|4-tc2|k3-k2|2 


(13) 


Notice that the coupling constant u(k;k 2 ;k 3 ) is zero under two conditions: (1) If all k allowed in '0(k) point at one 
single direction, say k; or, (2) if all k lie on one single ring, say |k| = k. For this reason, ordered states such as 
parallel rolls, hexagons, concentric rings, etc., do not have significant mean-flow couplings. Furthermore, the coupling 
constant t;(k;k 2 ;k 3 ) seems to have a pole at k 2 = ka. The real situation, however, is more subtle. Assume that 
ka = k 2 -I- q with q very small; then v ~ (e^ • k x q)(e^ • k 2 x q)(k 2 • q)/(c^( 7 ^ -I- q'^). For rigid-rigid boundaries 
(c^ = 10), there is no pole (since u ~ g) at g = 0. But for free-free boundaries (c^ = 0), a pole normally exists (since 
V ~ 1/g) unless q||k, or q||k 2 , or q _L k 2 . 

In this paper, we will mainly focus on two global quantities: One is the total convective current defined by 

J(t) = — f c?r^/>^(r,t) = ^ J(k, t) with J(k, t) = ^/;*(k, t)7/)(k, t); (14) 


the other is the total vorticity “current” defined by 

n{t) = jJ dru;lir,t) 

= Y /(ki+k2;k2)/(ki-f-k2;k3)^/>*(ki)^/>*(k2)V;(k3)V;(k4)4 i-l-k2,k3+k4- (15) 

ki,k2,k3,k4 

Notice that f(ki + k 2 ;k 2 ) ^ {kf — fc|)(ez • k 2 x ki). So n(t) = 0 if all the wavenumbers allowed in '!/'(k) point at 
one single direction ±k or lie on one single ring |k| = k. In other words, the vorticity current must be generated by 
couplings between modes of different k and k. From Eq. (|^), it is easy to derive that 

ro9tJ(k, t) = 2r(e;fc)J(k,t) 

- Y [ff(k-k2)+u(k;k2;k3) 

k2,k3 


'0*(k)^*(k2)'!/i(k3)^(k -I- k 2 - ka) -|- c.c. 


(16) 


In principle, this is the equation determining the structure of the convective current J(k, t) which, however, is beyond 
our present goal. Now applying the relations '(/'*(k) = 'i/’(—k) and u(k;k 2 ;k 3 ) = —u(k;k 3 ;k 2 ), and, exchanging the 
summation indices k —> —k, k 2,3 —> —k 2,3 for the g terms and k ^ k + ka — ka, ka ka for the v terms, one obtains 
from the above equation and the definition of J(t) that 

iToatJ(t) = ^r(e;fc)J(k,t) - Y 5 (ki • k2)Vi*(ki)V)*(k2)V’(k3)^/i(ki-h ka - ka). (17) 

k ki,k2,k3 
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This is the equation determining the total convective current J{t). Notice that the v terms vanish from this equation, 
which can also be derived directly from Eq. (|^) by converting the corresponding integral in Eq. (0) into a surface 
term. In general, the v couplings affect J(t) implicitly by modifying its structure J(k, t) unless, of course, u = 0. 

For stationary states, the convective current and the vorticity current are time-independent. This, however, is 
no long true if the state is spatiotemporal chaotic. For a spatiotemporal chaotic state, these two currents normally 
fluctuate in time around some well-defined averaged values: see Refs. @0- While the fluctuations appear chaotic 
in time, they are relatively small in comparison with their averaged values. For simplicity, we only consider the two 
corresponding time-averaged currents in our theory. We now introduce the time-average operator T defined by 

TF{t) =W)= ^lim ^ r dtF{t). (18) 

T^+oo 1 Jq 

Applying T to Eq. © yields 

y]r(e;fc)J(k,f) - (?(ki • k2)^*(ki)?^*(k2)?^(k3)V;(k4) ^ki+k2,k3+k4 0- (19) 

k ki,k2,k3,k4 

In the next several sections, we show how, under various assumptions, to calculate the time-averaged convective 
current of STC from this equation. The time-averaged vorticity current can be obtained with T acting on Eq. (HI)- 

For simplicity, we denote from now on J(k) = J(k, f), J = J{t) and fl = 11 (f). 

Finally we introduce the time-averaged structure factor defined by 

5(k) = J(k)/J with y]5(k) = l, (20) 

k 


and the corresponding averages 


(F)k = 05(k)F(k). (21) 

k 

With this notation, the first term in Eq. ( |T^ can be rewritten as (r(e))kJ. If the ^-dependence and the angular 
dependence in 5(k) can be separated, then it is more convenient to define 

nCO n27r 

S(k) = (2TT)^A~^S{k)^{a) with / dkkS{k) = l and / da^{a) = l, (22) 

Jo Jo 

where a is the angle between k and some reference direction. Here the discrete k lattice has been converted into 
a continuous one. So a proper phase factor has been taken into account. Notice also that d)(7r -|- a) = $( 0 :) since 
;§(—k) = 5(k). Now the corresponding averages with respect to S{k) and $( 0 :) are defined as 


poo p2'K 

{F)k= / dkkS{k)F{k) and {F)a= / da^{a)F{a). 

Jo Jo 


For F(k) = F{k, a), it is easy to see that (F)k = {F)k,a if a separation like Eq. (p^ holds. 


(23) 


III. A SIMPLE THEORY FOR PT 

In order to calculate the time-averaged convective current and the time-averaged vorticity current from Eqs. Jl^ ) 
and it is obvious that more information about the corresponding state is needed. We now present a simple model 
for PT, following a conjecture by Busse and his coauthors |Q, and calculate the corresponding convective current. 
For this task, we notice first that a PT state consists of many superimposed rolls with different orientations 
among which no particular direction is preferred. On the other hand, since this PT occurs immediately above the 
onset |5|Jl5|-|r7|| , the amplitude of the wave number k lies in the vicinity of kc- Our numerical simulations also 
revealed that the time-averaged structure factor is isotropic azimuthally. 

As a first attempt, we assume that a PT state is composed of many parallel rolls, whose wave numbers lie on a ring 
and whose amplitudes are equal. More precisely, we assume that 
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(24) 


M 


^/>(k, t) = tpo{t) 5k,q 


i=l 


where tpoit) and 4>{(3i,t) are the amplitude and the phases of those selected modes, ( 7 (~ k^) and /3i{t) are the amplitude 
and the angles of their corresponding wave numbers, while M is the total number of those modes, and a is the angle 
between k and some reference direction. From Eq. one finds that J{t) = = 54tpQ{t). One has 

M = 2, ipo = = const., and /3j = in with i = 0 1 for parallel rolls, and M = 6, ipQ = \fjJ%, 

t) = const., and Pi = in/3 with i = 0,1, • • •, 5 for hexagons |^. For PT, we take M +oo. Since J(t) in 
PT has a well-defined time-averaged value with small fluctuations p^ , we expect the same behavior for But 

we speculate that both the phases 4i{Pi,t) and the angular distribution {/3i(t)} are irregular in time, which leads to 
the spatiotemporal chaotic behavior in PT. From '!/;*(k, t) = ^(—k, t), one finds that if Pi{t) is selected, so must be 
n + Pi{t) with (j){n + Pi,t) = -(j){Pi,t). 

We now use this model to calculate the time-averaged convective current from Eq. di). While the first term can be 
easily evaluated as r(e, g) J, the second term is much more complicated. Notice that the condition ki -f k 2 = ka -|- 
imposes a very strong constraint on the available wave numbers on a ring. There are only three possibilities for the 
condition to be satisfied, see Fig. ||: (a) If ki -t- k 2 = 0, then ka -|- k 4 = 0; (b) if ki 7 ^ k 2 and ki -f k 2 7 ^ 0, then either 
ka = ki and k 4 = k 2 or ka = k 2 and k 4 = ki; or, (c) if ki = k 2 , then ka = k 4 = ki. It is more convenient to express 
these constraints in terms of their angles, which can be summarized, correspondingly, as: (a) 5a^^+'K5ai,a3+'K, (b) 
(l-^a 2 .ai-i- 7 r-< 5 a 2 ,ai)(^a 3 .ai< 5 a 4 ,a 2 +< 5 a 3 .a 2 ^a 4 .ai), and (c) ^ 02 .ai^a 3 ,ai< 5 ^ 4 .ai• Now inserting Eq. into the second 
term of Eq. (19) and applying -|- tt) = —(j){Pi,t) and these constraints, after some algebra, one finds that the 


second term is simply with 


9 1 2 J _ 

9M = 5(-l)(l - + 9[cosiP^{t) - Pj(t))]. (25) 

*,t=i 


Here V'o(^) and g[cos(/3i(t) — Pj{t))] have been decoupled, which seems reasonable. Since J = if one neglects 

the fluctuations of poit), one has then ~ J^. From Eq. (p^, this leads to the solution for the total 

time-averaged convective current 


J-r{e,q)/9M, (26) 

in addition to the conduction solution J = 0. This solution reproduces the known results for both parallel rolls 
with 32 = gi-l) + yil) and hexagons with ge = \[Qg{-\) + Ag{-\) + Ag{\) + g{\)\. 

For PT, we expect that Pi{t) and Pj{t) are uncorrelated in time and that fluctuations of Pi{t) can be neglected such 
that P/{t) ~ Pi{t) for any positive integer n. So we may make another approximation such that g[cos{Pi{t) — Pj{t))] ~ 
g[cos{Pi{t) — Pj{t))], where {Pi{t)} should be uniformly distributed between [0, 27r]. From this approximation and 
YYi {M/2n) /p ^ da as M ^ -l-oo, one finds that 

2 r 

goo=g{-l) + - dagicosa). (27) 

JO 

Using the explicit formula given in Ref. for free-free boundaries, one has finally that 

gpT = 0.855922 -h 0.0458144cr-i -h 0.0709326cr-^ (28) 

where a is the Prandtl number. Since typically g ~ fcc(l -I- qae), one finds that the time-averaged convective current 
J ~ r(e,g)/ 3 PT — ^/gpT in PT, recalling Eq. (||). For a = 0.5, this simple theory gives gpT — 1.2313. In comparison, 
we found gpp — 1.27±0.03 from our three-dimensional numerical calculations |l7|| . Considering all the approximations 
we have made, such a good agreement is very encouraging. 

However, this simple model apparently misses two important features of PT. The first is the lack of the vorticity 
current. Since all k lie on a single ring in our model, the vorticity current is identically zero from our discussions in 
the previous section. This, however, is not born out by our numerical calculations. Secondly, the structure factor 

from our numerical calculations has a finite width near its peak position, which leads to a significant reduction on the 

value of J as shown in Sec. So it is inaccurate to assume |k| = q for all k. To improve it, we now discuss a more 
physically sound theory. 
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IV. RANDOM PHASE APPROXIMATION FOR STC 


A. Convective Current 

We now consider a different approach to STC (both PT and SDC) in RBC. For this purpose, we notice that 
although the instantaneous patterns in STC are irregular and random in space MU, some spatial symmetries 
can be restored if the system is averaged over a very long tim e. For example, for a laterally infinite system, it seems 
natural to assume that the time-averaged quantity is uniform in space and that the two-point correlation 

function, defined as C'(ri,r 2 ) = '0(ri, i)'!/'(r 2 , f)/'0^(ri, f), is translation invariant, i.e., C(ri,r 2 ) = C{vx — r^)- Then, 
one finds from Eq. ( p^ that t) = J, and, from Eqs. and d^ ) that 

= J(5ki,k2'5(ki). (29) 


It is also easy to check that 5(k) is just the Eourier component of C'(r). 

To understand the physical implications of the above result, we assume that 

'0(k,t) = /5(k,t)e*'^(‘"’*\ 


(30) 


with both /5(k, t) and (^^(k, f) real. Since = ■!/i(—k, t), one has /5(—k, t) = p(k, t), 0(—k, t) = —0(k, t) and 

J(k) = While both experimental measurements and numerical calculations suggest that the amplitude has 

a well-defined time-averaged value with small fluctuations, the phase seems to have a rather complicated, irregular 
time-dependence. If the amplitude and the phase can be assumed to be uncorrelated in time, the above result indicates 
that phases of different modes are totally uncorrelated in time. Since the Jki k 2 factor in Eq. ( p^ can be represented 
by 


2 A-1 


iJki,k2 = (27r) A 


1 


(5(ki - k 2 ) = lim (27r^^/A) exp[--(ki - k 2 ) 


(31) 


we hence adapt a random phase approximation (RPA) to STC in which 

exp[-4(ki,t) -hi(^(k 2 ,t)] = exp[-i(ki -k 2 )^^^]. (32) 

Here is a length determining the correlation between phases of different modes. We expect that = [A/27r]^/^ ^ 

-|-oo for a laterally infinite system but, as we will show, this limit should be taken only later on. 

Now, since the phases are random in time, one may further approximate a four-phase correlation by products of 
two-phase correlations, i.e., 


exp[-z^(ki,t) - i(j){k2, t) + icjiikzA) + f^(k4,f)] 

~ exp[—i(/)(ki, t) — i0(k2, t)] exp[i(j){k 3 , t) + #(k4, i)] 
-|-exp[-i(^(ki,t) +i^(k 3 ,t)] exp[-z(^(k2,f) +i^{ki,t)] 

+exp[-i^{ki,t) Ai^ikiA)] exp[-i(^(k2,t) +i${k 3 ,t)] 

2f2 


= exp 


~ ^(ki + k2)^^0 — -(ka -I- k4)^^^ 


exp 


-i(ki-k3)2e^-i(k2-k4)2C^ 


+ exp 


-i(ki-k4)2^|-i(k2-k3)^^^ 


(33) 


Applying this and neglecting the correlations of J(k,f) such that J(ki, t) J(k 2 ,f) ~ J(ki)J(k 2 ), one may rewrite Eq. 
(p^) in a continuous Eourier space as 


(r(e, fc))kJ = 




c?kidk2dk3C?k4 


(27ry 


A 

\2c2 1 


(5(ki +k2-k3- k4) 5(ki • k2) 




( 34 ) 
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Since +oo, we may take g(ki • k 2 ) = 5 (—1) for the first term on the right hand side. Then, after some algebra, 

one finds that 


_ {rie,k))k _ 2 (r(e,fc))fc 

^ 9csd !Joo 


in addition to the conduction solution J = 0. Here A = 27r^^ has been used, ^(k) is assumed to be azimuthally 
uniform, and g^o is defined in Eq. (^). Notice that this result is valid for both PT and SDC, but the exact values of 
5 oo and (r(e, fc))fe depends on the boundary conditions. While the value of g^o for PT has been given in Eq. (28), the 
corresponding value for SDC is 


= 1.1319 + 0.0483cr-i + 0.0710cr" 


(36) 


which is obtained by interpolating the data for g{cosa) given in Ref. |^| for rigid-rigid boundaries and by integrating 
the consequent fitting function. 

It is worthwhile to point out that if one takes the limit -l-oo as early as in Eq. (33), one misses the factor 2 in 
Eq. d^). To understand this, one should notice that Eq. ( |^ ) reduces to ( 5 ki,-k 2 < 5 k 3 ,-k 4 + iJki,k 3 < 5 k 2 ,k 4 + iJki,k 4 < 5 k 2 ,k 3 
under the limit -too so that the constraint ki -|-k 2 = ka -|-k 4 in Eq. ( |l^ is automatically satisfied. Consequently, 
the delta function (27r)^44“^(5(ki-|-k2 —ka —k 4 ) in Eq. ( ^ ) is not necessary and can be replaced by 1. Considering that 
this constraint is intrinsic in our problem and physically needed, we believe that such a replacement is not justified. 
In our calculation, the limit > -Poo is taken only at the last step of the calculation, which makes the constraint an 
important feature in Eq. (^^. Our approach is also justified by our numerical results. It is obvious from Table | and 
H that this factor of 2 improves the agreement between our theory and our numerical results. 


B. Vorticity Current 


We now use the RPA to calculate the time-averaged vorticity current. From Eqs. (p^ and (133), one finds that 


n = 


(^ 


dki dk2 dks dk4 / (ki -P k2; k2) / (ki -P k2; ka 


(27r) 

A 


A(ki -Pk2 -ka -k4) 




(37) 


where, from Eq. 


/(ki -P k2; k2)/(ki -P ka; k3)(5(ki -P ka - ka - k4) 

(fci - - fc|)(ez • ki X k 2 )(e^ ■ ka x k4) , u _ u 

4a2 [|ki+k2P+c2][|ka+k4|2+c2] 1+ 2 3 


k4). 


(38) 


Notice that there is a singular point ki -P k 2 = 0 in the above expression for free-free boundaries: Since = 0 from 
Eq. (|^, the value of the above expression depends on how the point is approached. On the contrary, it is smooth 
everywhere for rigid-rigid boundaries, in which case (? = 10 from Eq. (|^). 

The evaluation of D is rather complicated, which we present with some details here. We assume that .S(k) is 
azimuthally uniform, i.e., J(k) = 2TrA~^ JS{k) from Eqs. ( po|) - Clearly D consists of three terms. With 

q = ki -P k 2 , the first term, after some algebra, can be reduced to 




_/ da — _ 

(2^)4 4^2 J + c2)2 


dki S'(fci)( 2 q • ki - q'^){e^ • ki x q) 


X 


dks S'(fc 3 )( 2 q • k 3 - q'^){e^ ■ ks x q). 


(39) 


It is easy to see that the angular integrals over ki and ks are both zero, so Di = 0. By interchanging k 3 ^ k 4 , one 
can show that the second and the third term are identical, so D = 2 D 2 = 2 D 3 . With q = ki — ks, one finds that 














(27r)4 20-2 


n = 


dkidk2S'(/fci)S'(/c2) 


(/c?-fc|)(e^-ki X k2) 


[|ki +k2|2+c2]2 

X J dqe~‘^ - A:| - 2 q • (ki + k2)] [e^, • ki x k2 - • q x (ki + k2)] 

p r r 

ym^ 


4tT2 


dkikiS{ki) / dk2k2Sik2)A{ki-k2;c^), 


where we have used A = 2Tr^'^ and 


(40) 


A{ki; k2-, c^) 


{kl-klfklkl 


-^{kl-klf 


da 


sin^ a 


TT [kf + A:| + 2 kik 2 cos a + P]'^ 


fc2 + k. 


\/{kf - fc|)2 + 2P{kl + fc|) + 


- 1 


(41) 


For free-free boundaries, since P = 0, one has that 

A(fci;fc2;0) = _ fc2| [^2 + _ 1^2 _ ^2|]^ (42) 

which has a second-order singularity at ki = k 2 and is due to the singularity in Eq. (|^). In comparison, the function 
A(fci; k 2 ',P) is analytic everywhere for rigid-rigid boundaries with P = 10. While A(fci; k 2 \ 0) ^ for free-free 

boundaries, correspondingly A(A:i; ^ 2 ; P) (kf — fc|)^ for rigid-rigid boundaries. As we will show in the next section, 
this has a significant consequence to the properties of STC. 


V. SCALING RELATIONS IN STC 


A. General 

To evaluate the convective current J and the vorticity current O, one must know the structure factor S(k) which, 
however, is beyond our present theory. We thus turn to phenomenological arguments. We define a two-point correlation 
length as 

i2=[{k^)k-{k)l]~^'\ (43) 

Then we assume that the structure factor satisfies the following scaling form 

kS{k) = ^2mk - k^aP^2], (44) 

where k^ax is the peak position of kS{k) and P{x) is the scaling function satisfying J^^dxJ-(x) = 1. [Since fc > 0 
in kS{k), the lower limit for J^{x) is —kmax^ 2 , which we approximate by — 00 .] Inserting k = kmax + x^ 2 ^ and Eq. 
( ^ into Eq. (^), one gets that {x'^)x — {x)^ = 1, where we have used the notation 

/ CX> 

dx T{x)F{x). (45) 

-00 

It is also easy to see that {k)k = kmax + ^2 ^ ^ ■ 

For very large ^ 2 , one may take k = kmax + x^ 2 ^ in Eqs. ( ^ and ( ^^ and expand the results in order of 1/^2- It 
is easy to find from Eq. (||) that 

{r{e,k))k = r{e,kmax) - b{p - l){x}xkaeja - - l){x^)xeo/^2 + 0{l/e2), (46) 

where b = kmax/kc- We expect that 6 ~ 1 -|- 6ie in STC. So for small enough e, one has that 

{rie,k)),^e-{x^)a^eo/e2, (47) 

and, from Eq. (^^, that 
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J' 


9o 


e-(x“)4 

S2 


21 


(48) 


with goo = gpT or goo = gsDC depending on the boundary conditions. This expression depends on two unknown but 
experimentally measurable quantities ^2 and {x'^)x = 1 + If J^{x) is symmetric, then {x)x — 0 and {x'^)x = 1- 
In general, however, one has {x‘^)x > 1- Notice also that although mean-flow couplings are not explicitly present 
in Eq. (H), they affect the value of the convective current via the structure factor 5(k) [see Eq. ®] and the 
two-point correlation length ^ 2 . For PT, since ^2 — 0. the contribution from ^2 reduces the value of J quite 

significantly. This feature is absent in the simple theory presented in Sec. III. 

One may evaluate the vorticity current in the same way. From Eqs. (^) and ( ^ ) , one gets for free-free boundaries 
that 


n : 


ym^c ^ 
4cr2 ^2 


([xi - a;2|)a 


(49) 


where we have used kmax ~ kc- The quantity (jxi — X 2 \)xi,x 2 is related to the width of the scaling function T{x) and, 
since {\xi - X2\)x^,x■2 < \/{{x\ - X2Y')x^ we expect (|a;i — X2\)xi,x2 to be of order of unity. Similarly, we find 

from Eqs. (^^ and © for rigid-rigid boundaries that 




„2 u2 t2 

ym^c ^ 


2k‘i + (? 




- 1 


2a2 el 

Clearly, by phenomenological arguments, we can express J and Q in STC in terms of measurable quantities. 


(50) 


B. PT 

For PT, we further assume power law behaviors for the two-point correlation length, the convective current and 
the vorticity current such as 


Then, from Eq. 


Recalling Eq. 


C 2 ~ ^ 2,06 J ~ and 
]), we find the following scaling relation 

A = 2/i 12 . 

one obtains that 


J ■. 


2 

gpT 


S2,0 




Since J is positive by definition, the values of the exponents satisfy 

p, = l, j^>l/2 and \ = 2 + v>b/2 


in PT. 


It is very likely that v = 1/2, hence, A = 5/2. If so, then one finds from Eqs. (48) and ([49|) that 


Jo = 


gpT 




2,0 


and = 

4cr ?2,0 


(51) 


(52) 


(53) 


(54) 


(55) 


which depend on three phenomenological parameters ^ 2 , 0 j {x'^)x and (jxi — X 2 \)xi,x 2 - If *2 > 1/2, then Jq = 2/gpT 
since the term in Eq. (^) contributes only to the leading correction to scaling. It is interesting to notice that the 
amplitude equations coupled with mean-flow predicts for free-free boundaries that Q ~ for almost perfect 
parallel rolls. 

We now verify our predictions for the power laws in PT by our numerical solutions. We have carried out large-scale 
numerical calculations of the three-dimensional Boussinesq equations under free-free boundaries for fluids of ct = 0.5 

f j. We have confirmed in Ref. 
one can see that our theoretica 


17| ] that the structure factor in PT satisfies the scaling form (^^. From Table 
and our numerical results are in very good agreement for the exponents. The 
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scaling relation Eq. ( p2D is confirmed within our numerical uncertainties. The comparison between the corresponding 
amplitudes, however, is only moderately successful. Calculations of ^ 2 , 0 ; {x^)x and {\xi — X 2 \)xi,x 2 are obviously beyond 
the present theory, so we take our numerical result for ^ 2 , 0 - Since our numerical results for {x^)x and {\xi — X 2 \)xi,x 2 
are too sensitive to the large value cutoff to be me aningful, see disc ussions in the next section, we assume equalities 
in {x^)x = 1 + {x)l > 1 and (|a:i - X 2 \)xi,x 2 < \/((a^i - X 2 y)xi,x 2 = From Eqs. (||), (|) and @, one gets 
Jo ~ 0.972, which is about 20% larger than the numerical value. A non-zero value of (^x will apparently reduce 
the theoretical value of Jq in the right direction. It is worthwhile to point out that, since ^ 2,0 — (3/2)^o |17|, the 
value of Jo is reduced significantly owing to the finite width of the power spectrum. On the other hand, one finds 
that Oo ~ 454.7(|a;i — X 2 \)xi,x 2 = 643.0 as an upper bound, which is about ten times larger than our numerical 
result. Nevertheless, we note that while our theory is based on the two-dimensional GSH equations, our numerical 
calculations are done for the three-dimensional Boussinesq equations. Although the former is very good in reproducing 
qualitative features of RBC, it may not be quantitatively accurate in modeling RBC |^,|^. So one should be cautious 
in comparing the results from the GSH equations with those from real experiments or those from numerical calculations 
with hydrodynamical equations. 


C. SDC 

The situation for SDC, however, is more subtle since the roll-to-SDC transition occurs at a positive temperature 
€t P-p^. Consequently, several competing scaling scenarios are possible in SDC, the choice of which depends on the 
character of the transition. A more thorough examination on the issue will be presented elsewhere [p^ . We mention 
that the same power laws as Eq. can be defined for SDC. But instead of Eq. (p^, one finds from Eqs. ( |^ and 
(pl|) the following scaling relation 


A — 2^ -\- 2z/. 

By the same arguments leading to Eq. (^J), one gets that 


(56) 


Ai = l, J^>l/2 and \ = 2 + 2v>2, in SDC. (57) 

The conclusion that different scaling relations hold for PT and SDC can be traced back to Eqs. (|^) and (p^ ) via the 
different behaviors of A(fci; ^ 2 ; c^) at fci = ^ 2 - 

In order to test our theory of SDC, we have carried out systematic numerical studies of SDC with the GSH equations 
p^ . For simplicity, we take gicosa) = g as a constant so, from Eq. (p^, = ^g- For numerical convenience, 

following Refs. [||Q, we rescale the GSH equations such as 


^r', t (4To/fc2^^)t', '0 ^ {kc^o/‘2^/g)ip', 

e - a - (eoV4To)a', ^ 


C ^ {9mTokl/g)C,', 

9m {9ill^Tlkl)g'^, 


which leads to the rescaled GSH equations 


(58) 


dt.^' + g'^\J' = e'-(V'' + l)2 


fi 




V'^C' = e. 


V'(V'>') X VV' 


(59) 

(60) 


where U'(r') = V'C'(r',t') x e^. Now the time-averaged convective current 
current (p® are rescaled into 


and, 


J( 


SDC 


/ 

e — 


4 (a 


cD 




SDC 


2 ct '2 


2-he' 


n 


V4c'^ -h d‘ 


- 1 




and the time-averaged vorticity 


( 61 ) 


(62) 


From Eq. (||), one finds that e = 0.3594e' for rigid-rigid boundaries. In principle, for a given cr and a suitably chosen 
g, the parameters g'^, a' and c'^ are determined by Eqs. (||) and (|^). Again following Ref. |^, we simply choose 
g'^ = 50, a' = 1.0 and c'^ = 2.0. Details of the numerical studies of SDC are presented elsewhere p^. 
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One crucial assumption in our theory of SDC is that the structure factor S{k) has a scaling form like Eq. (|4^). So it 
is very important to verify this assumption. In the insert of Fig. the results for k'S'{k') for e' = 0.55, 0.65 and 0.8, 
corresponding to SDC states, are plotted. [The structure factor S'{k') is nomalized by dk' k'S'{k') = 1.] To check 
whether a scaling form like Eq. ( [44| ) holds, we take the two-point correlation length ^2 from our numerical results and 
choose to give the best fit to scaling. For each e' within 0.55 < e' < 0.8, we hence find a corresponding function 
lF(x) of SDC, which is shown in Fig. H. As one can see, all the data collapse into one single curve. The scattering of 
the data near is due to our numerical uncertainties and is within the corresponding error bars. So the existence 
of a scaling form of kS{k) is verified within our numerical uncertainties for SDC. 

We now compare our numerical results for J' and D' with those from our theory. Theoretical results are presented 
in Eqs. and (|6^. We fit our numerical data with power laws such as ^2 = ~ J' = ~ ^'cY ~ 

and Vt' = ng(e' —e[,)^ with e[, = 0.002, see Ref. |^. The non-zero value of e[, is likely due to finite-size effects. In Table 
H, we summarize both theoretical and numerical results for J', D' and for SDC. We actually put /i = 1 in our fitting 
of J', so the agreement with this is trivial. The inequality for the theoretical value of J'^ is from {x^x > 1- Since the 
calculations of ^2 0 ^ beyond our theory, we use the corresponding numerical results in calculating fig and A. 

Clearly the scaling relation Eq. (^) is approximately verified. The prediction for J' is very good. The prediction 
for the value of Dg, however, is larger than the corresponding numerical result by a few magnitudes. The cause for 
such a big discrepancy, at present, is not clear to us. Considering that ujz{v,t) has a highly localized structure in 
real space it is possible that our numerical calculation is not long enough to sample all the phase space. It is 

also possible that assumptions in our theory are not sufficient to describe the behavior of fl. In comparison with the 
situation in PT, which is discussed in Sec. 0(B), the success of our theory in describing SDC is not as satisfactory. 
Further improvement of it is obviously valuable. 


VI. IS POROD’S LAW VALID? 


In phase ordering, a sharp interface exists between domains of different phases. Consequently, the real-space 
correlation function C(r) is proportional to rjL at short distances, where L is a characteristic length of the system 
pSf . Then the corresponding structure factor, which is the Fourier transformation of C'(r), behaves like S{k) ~ 1/Lk^ 
for large k in two-dimensional space. This large k behavior of S{k) is known as Porod’s law p7| , p8|| . It is easy to check 
that the two-point correlation length defined in Eq. (^) is very sensitive to the large k cutoff if Porod’s law is valid. 
As a result, other criteria are needed to define a better behaved characteristic length, say L, of the system. 

For the convective patterns in RBC, smooth interfaces are always present between hot, rising fluid and cold, sinking 
fluid. In the ordered states, the patterns can be described by a few sine or cosine modes. Correspondingly, the 
structure factor consists of only several sharp peaks. Porod’s law is not relevant in this case. But in STC, an infinite 
number of modes are excited, including those large k modes. Then, a natural question can be raised: Is Porod’s law 
valid in STC? Considering that the shape of the interface between different domains appears random and the motion 
of it seems chaotic, an intuitive argument is rather difficult. In this section, we present our efforts in this direction. 

To start, we take the scaling form Eq. (|4^ ) of kS{k) but replace ^2 with a characteristic length in case ^2 is 
cut-off dependent. If ^2 is well-defined, from Eqs. (|4^) and (Q (with ^2 replaced by ^s), one gets that 


« 6 \/YYx - {x)x- 

So and ^2 are identical up to an overall constant. Recall the following formulas 

CX> 

gzk r _ _|_ 2 ^ i^Jmikr) cosma, 


(63) 


(64) 


Jo(x + y) = Jo(x) Jo(y) + 2 ^ (-l)™Jm(x) Jm(y), (65) 

m—1 

where a is the angle between k and r, and Jm{x) is the m-th order Bessel function. It is straightforward to show that 

P 00 

C{r) = J dkkS{k)MkT) = Mk^axr)Co{.r/£,s) + 2 ^ {-IT Jm{k^axr)C^{r / £,s), (66) 


where 
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^+oo 


dxT{x)Jm{xr/^s)- 


(67) 


Since Jm{y) ~ y™, one has that Cm{y)ly^ ^ 0 as y ^ 0 + for alH = 0 , 1 , • • •, to — 1 . So, neglecting the possible 
presence of singularity, we assume the following expansions 

OO 

Cm{y) = y™ X! y" y > 0- (68) 

i=0 

[One cannot apply the small y expansion of Jm{y) in Eq. ( | 6 ^ since, for any fixed r/^s, the integral is dominated by 
those x’s such that xr/^s is not small.] Since Jo(y) ^ 1 as y ^ 0+, one finds that Coo = dxT{x) = 1. Now it 
is easy to see from Eq. (|^ that 


C(r) « 1 + Cqit/^s for kmaxX ^ kcT < 1. 


(69) 


While the constant term contributes an unmeasurable i5(k) to 5(k), the linear term leads to the Porod’s law, i.e.. 


S'(fc) - for 


(70) 


It is worthwhile to mention that the 1/^s dependence is as important as the l/fc^ dependence In phase 

ordering, a large k cutoff exists so that Porod’s law is valid for those fc’s smaller than this cutoff |2^ ]. It is not clear 
whether such a large cutoff exists in STC. One possibility is that this cutoff exists and is of the same order of kc, in 
which case Porod’s law is limited to such a narrow range in k space that verification of it is almost impossible. 

Assuming ^2 is well-defined, we plot ^ 2 k^S{k) vs. k for both PT and SDC in Fig. The data for PT are obtained 
from our numerical solutions of the three-dimensional Boussinesq equations jp^ , evaluated at the mid-plane. The data 
for SDC are from our numerical calculations of the GSH model [^. As one can see, the value of ^ 2 k^S{k) in PT seems 
to approach a constant at large k, insensitive to the exact value of e. So Porod’s law might be valid in PT. But the 
value of ^ 2 k^S(k) seems to increase for large k in SDC! However, it is known the GSH model introduces an artificial 
short-ranged (hence large k) cross-roll instability so the large k behavior in the GSH model might be different 
from those in real systems. Furthermore, owing to the finite grid size used in our numerics, we are not sure how 
numerical noise might affect the large k behavior in both PT and SDC. For this reason, we believe that more accurate 
data are needed for a definite conclusion. Even so, one sees immediately how sensitive the two-point correlation length 
(2 defined by Eq. (|^ could be to the large k cutoff. So it is useful to define a less sensitive characteristic length, say 
^S, of the system. One obvious choice is the inverse of the full width at the half peak (FWHP) of kS{k). Since one 
can easily find a function iF{x) satisfying Eq. ( ^ ) for each ^5 (replacing ^2 with ^ 5 ), this provides the easiest way 
to check whether a scaling form exists. If the system is inside the scaling range, all JF(x)’s so defined should collapse 
into a single curve. One must, of course, normalize S{k) by f dkkS{k) = 1 first. But this normalization is much less 
sensitive to the large k cutoff than ^2 is. As shown in Eq. (^), if ^2 is well defined, ^5 and ^2 are simply proportional 
to each other inside the scaling range. This is not true if the system is outside the scaling range. 


VII. DISCUSSION 


Our phenomenological theory for STC in RBC depends on two basic assumptions. In Sec. IV, we assume that 


the time-averaged two-point correlation function is translation invariant in real space and we hence adapt a random 
phase approximation to STC. In Sec. ^ we further assume that the structure factor satisfies a scaling form such 
as kS{k) = ^ 2 iF[{k — kmax)i 2 \- In comparison with similar scaling forms in critical phenomena, critical dynamics 
and phase ordering we find it necessary to replace k with k — kmax in the scaling form. The physical origin 

of this replacement is due to the fact that patterns in RBC have an intrinsic wavenumber, which is close to kc- 
By the same reason, we find it necessary to seek the scaling form of kS{k) instead of S{k), where the k factor 
comes from dk = kdk da in two-dimensional fc-space. The existence of the scaling forms in critical phenomena and 
critical dynamics is rooted in the scaling invariance of long wavelength fluctuations in the system and is associated, 
respectively, with a fixed point in renormalization group theory [Q. Its physical origin in STC is yet unknown. In 
Sec. ^ we have confirmed the scaling form of S{k) within our numerical accuracy. Since fc > 0 in kS{k), the lower 
limit for the scaling function iF{x) is —kmax^ 2 , which is e dependent. So the violation of scaling is almost certain for 
very small k. We cannot rule out from our numerical data that this scaling form might also be violated for very large 
k. It is not clear currently in what range the scaling form is valid. 
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As we discussed in Sec. Vi, the two-point correlation length ^2 is cutoff dependent if Porod’s law is valid for 
STC in RBC. In principle, there is another disadvantage to choose ^2 as a characteristic length. It is easy to see 
from Eq. that {k)k = kmax + so {k)k is shifted from k^ax by {x)xl^ 2 - Because of this, an unknown 

parameter {^)x is introduced in Eq. (|48|). This {x^)x parameter can be easily removed by defining a new length 
Cs = [{{k — kmax)'^)k]~^^'^, instead of Eq. (^. Then one simply has {x'^)x = 1 if .^2 is replaced by in Eq- @)- 
In practice, however, our numerical data are not accurate enough to determine kmax precisely. Consequently, there is 
no practical advantage for us to use instead of ^ 2 . This may not be true for experimentalists since their data are 
much more accurate. Of course, it is also to be tested whether the structure factor can satisfy a scaling form like Eq. 
with respect to so defined. 

In summary, we present a phenomenological theory for STC in RBC. We calculate analytically the time-averaged 
convective current J and the time-averaged vorticity current Q in both PT and SDC as functions of e and ^ 2 - Our 
theory is successful for both PT and SDC, despite the need for a better quantitative result for fl in SDC. We believe 
that our theoretical results will be useful in understanding the complicated behavior of STC in RBC. We also believe 
that our theory provides a new approach to STC and also raises some interesting questions. For example, how 
can one calculate the structure factor S{k) and the two-point correlation length ^2 analytically? Is it possible that 
certain global quantities in STC form a complete set in the same way as temperature, pressure and density do for 
thermodynamic systems? Can we derive some effective variational principle in terms of global quantities? How far 
can we apply the ideas in critical phenomena to study STC? Since our assumptions are quite general, it will also be 
interesting to see whether our theory can be generalized to STC in other systems 
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TABLE 1. Time-averaged convective current J « Joe^, time-averaged vorticity current H and two-point correlation 

length ^2 ~ ^ 2 ,oe~'' in PT with a = 0.5. For theoretical result of iz, we assume equality in Eq. (54). See also discussions in Sec. 
V(B). 



T 

V 

A 

^2,0 

Jo 

Ho 

Numerics 

1.034 ± 0.025 

0.472 ± 0.016 

2.55 ±0.10 

0.82 ± 0.04 

0.787 ±0.019 

70.1 ± 1.0 

Theory 

1 

1/2 

5/2 

— 

0.972 

643.0 


TABLE H. Time-averaged convective current J' « time-averaged vorticity current Q,' « and two-point 

correlation length ^2 ~ ^ 2 ,oc'~'' in SDG, with g'm ~ 50, a' — 1.0 and c'^ = 2.0. For numerical results, we actually use e — e(. 
with e[, = 0.002 instead of t for data fittings. For more details, see Ref. [26]. 



h- 

V 

A 

^2,0 

Jo 

4 

n'o 

Numerics 

1 

0.72 ±0.05 

3.0 ±0.1 

6.8 ±0.2 

0.64 ±0.02 

2.9 ±0.9 

(3.0 ±0.2) X 10"* 

Theory 

1 

> 1/2 

3.4 ±0.1 

— 

2/3 

> 8/3 

3.7 X 10“^ 
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FIGURE CAPTIONS 


Figure |^. Allowed configurations of wavenumbers satisfying ki = k 2 = ks = ki and ki + k 2 = ks + k 4 : (a)ki + k 2 = 
ka + k 4 = 0 ; (b)ki ^ k 2 and ki + k 2 7 ^ 0 ; or (c)ki = k 2 = ka = k 4 . 

Figure |[ A plot of k'S'{k ')/^'2 vs. x = {k' — k '^^^)^'2 for 0.55 < e' < 0.8 in SDC, showing scaling and the scaling 
function !F{x) defined in the text. The scattering of the data is within our numerical uncertainties. Insert: The 
time-averaged function k'S'{k') vs. k' for e' = 0.55, 0.65 and 0.8 in SDC. 

Figure |[ Plots of ^ 2 k^S{k) vs. k in (a) PT and (b) SDC. The error bars are plotted only for (a) e = 0.05 in PT 
and (b) e' = 0.65 in SDC. 
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